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Abstract. Asymptotic forms of the Hilbert-Scmidt and Hilbert norms of 
positive definite Toplitz matrices Qm = (p(j — k))j j^_ Q as N — » oo are de- 
termined. Here b(j) are consequent trigonometric moments of a generating 
non-negative mesure dcr(9) on [— 7r,7r]. It is proven that cr(8) is continuous if 
and only if any of those contributions is o( N ) . Analogous criteria are given for 
positive integral operators with difference kernels. 

Obtained results are applied to processing of stationary random signals, in 
particular, neutron signals emitted by boiling water nuclear reactors. 

1. Introduction 

This paper is motivated by the studies [V(tM01| devoted to the problem of 
detection of hidden unstable components in random neutron signals measured in 
boiling water nuclear reactors. We assume that a signal from a monitored system 
forms, during a sufficiently long time interval, a real- valued stationary random 
process £ (t) , t E Z, with discrete time, such that the means 

<*(*)> =o, {e (<)>=i. 

The correlation function of such a process 

b (t) := (£ (t) £ (0)) = (t + t') i (<')) >t,t'eZ, 
is a real-valued sequence, which admits the representation 

(f) b(t)= [ exp (it6) da (6) , 

J — 7T 

where a (9) is a non-decreasing bounded function on [— w, ir] R67 . By virtue of 
our assumptions 

b{0) =<t(tt) -cr(-Tr) = f , 
and for any #i, 6*2 such that < 9\ < 02 < 7r, we have 
(2) a{e 2 )-u{e l )=a{-e l )-a{-6 2 ) . 
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In general, the spectral distribution function a (9) which determines the process 
correlation function by Q, can be split into a sum 

(3) a (9) = a c (9) + a d (9) 

of a continuous non-decreasing function a c (9) and a non-decreasing step function 
id (9), and such a representation is unique up to constant contributions in a c (9) and 
<Jd (9) |R67] . Notice that both functions, a c (9) and ad (9), satisfy the condition J3J). 
Actually, the problem, formulated in |VGM0l| . was to find, in a real-time operation 
mode, whether the spectral distribution function of the random signal a (9) contains 
or does not contain a non-trivial component ad (9). We will call here for brevity the 
random process (signal) £ (t) stable if its spectral distribution function a (9) is 
continuous and unstable otherwise. In other words, the process is unstable if and 
only if 

a d (tt) - a d (-7r) > . 

The main task of the present work is to formulate criteria of stability of the process 
in terms of its correlation function b (t) . As a main tool we use the sequence of 
positive definite Toplitz matrices Qn = (b(j — k))j k = Q . This paper is organized in 
the following way. 

In Section 2 we find the principal asymptotic contribution of the Hilbert-Schmidt 
norm ||Qat|| 2 with N — > oo, and show that the process is stable if and only if this 
contribution is o (A). 

Section 3 contains a similar criterion, but with the Hilbert norm ||<3./v|| instead 
of ||Qjv|i 2 - We prove here that if A — ► oo, then ||Qjv|| = m • A + o(N), where m is 
the maximal jump of a (9). 

In Section 4 both criteria are generalized for continuous time processes or for 
positive integral operators with difference kernels. 

In Section 5 the efficiency of the above stability criteria for signal processing 
is discussed and illustrated by application to simulated signals and real neutron 
signals emitted by the Forsmark 1&2 boiling water reactor. 

2. Asymptotic form of the Hilbert-Schmidt norm of truncated 
correlation matrices and the stability criterion 

Let us denote by {Qn} , A = 1, 2, . . ., the sequence of Toplitz matrices (b (j — fc))^Z^ 
and let ||<3jv|| 2 be the Hilbert-Schmidt norm of Qn- 

(N-l \ 5 / N-l \ 5 

J2b 2 (j-k)\ = A& 2 (0) + 2^(A-fc)fe 2 (A : ) . 
j,k=0 J \ k=l ) 



Our assumptions imply that 
b 2 (k)-- 



exp (ik9) da {9) 



2 

2 



< r (0) = l . 



Therefore \\Qn\\ 2 <N. 

Theorem 2.1. The process £ (t) is stable if and only if 

jiml||Qiv|| 2 = 0. 
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Proof. Introduce 



■ 2 1 N{6-6') 
N 2 sin 2 \ (0 - 0') 



and 



f sin 2 IN (9 - 6') 

*n(0) = / V -da d (9') = 

K ' J N 2 sin 2 \{6-e>) 

1 ( sin 2 iiV (0 - 9 a ) t sin 2 ±iV(0 + a ) 



2^ \N 2 $m 2 l(9-9 a ) N 2 sin 2 \{9 + 9 a ) ) ' 

where {±0 Q } is the set of jump points of a (9) or, what is the same, of the points 
of growth of a d (9) , 

^m a = a{9 a +0)-a (9 a -Q)=a d {-9 a + Q)-a d (-9 a - 0) . 
Due to and Q, we have 

rf-i 

j,fc=0 



7T 7T 



iV-1 



( 5 ) ^//[E ex P C C? - fc ) ( fl - ] ^ da (9) = 



v i,fc=o 



$jv (0) da c (0) + 2 / $ w (0) (0) + / * w (9) da d (9) 



Let 



w (6) = max ^ |cr c (0) - ct c (0')| 



The continuity of cr c (0) implies lo (5) 0. 
Since 



sin 2 iJV (0 - 0') 

(6) — — < 1 

{ } TV 2 sin 2 ± (0-0') " 

and 



2 7T 

| sin x | > — \x\ , |x| < — (mod27r) 

7T 2 



then, for any < S < -, we have 

sin 2 ±7V (6» - 0') 
TV 2 sin 2 ± (0 - 6»') 



$v (0) = " 2 Y"~ P; da c (0') 



|0-0'|<<5 (mod27r) 

1 at — ft' 

>-da c (0') < 



sin 2 \N (0 - 0') 
TV 2 sin 2 | (0 - 0') " 



0-0'|>(5 (mod2ir) 



7T 2 



(7) 2 ui(6) + — {a c (ir)-a c (-ir)}<2uj(S) + 



N 2 L L w yj " w 7V 2 <5 2 
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Therefore, by an appropriate choice of S and TV, the first two integrals on the right- 
hand part of (JjjJ can be made arbitrarily small. Hence, these integrals tend to zero 
as N — > oo. 

Rewrite now the third integral on the right-hand side of J5| in the form 

1^ / sin 2 \N {9 a - 9 a >) sin 2 ljV(6> a +6> a ,) \ 

2 q ^ q mQmQ ' ^ sin 2 \ (9 a -B a .) + sin 2 \ (9 a +6^) ) ' 

and assume that o~d (9) has a finite number of jumps, then we can conclude, taking 
into account JBJ), J7J) and the inequality 

m a = a d (vt) - a d (-n) < a (ir) - a (-%) = 1 

a 

that 

(8) ^^ = \^<- 

a 

Thus, 



(9) 



maxm Q < JlE m « = J im F IIQjvII 2 ^ 



^/maxm^yi 



X) m Q — /nraxTOo, 



The generalization of © and, hence, of @ for the case of a d (9) having infinitely 
many jumps can be obtained by continuity using the standard method employed 
in the next section. □ 

The condition of Theorem 12. II in special cases can be specified. 

Proposition 2.2. If the process is stable and its spectral distribution function a (9) 
satisfies the Holder condition: 

(10) \a (0i) - a (0 2 )| <A\9\ — 9 2 \ v ,0<v<l, - w < 6 X , 6 2 < tt, 

then 

\\Q N \\ 2 = 0(n&) . 

N— >oo \ / 

Proof. Using the estimate (J7J) and the assumptions of the proposition we deduce 
that 



\Q N \\ 2 < s j2AS"N^+^, 0<*<~. 
By the minimization of the latter estimate in 8 we obtain that 

\\Qn\\ 2 < Vin^ N^h . 

□ 
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3. Asymptotic form of the maximal eigenvalue of a truncated 

correlation matrix 

Let us denote by A m (N) the maximal eigenvalue of a positive (i.e., positive 
definite) matrix Qn, A m (N) = ||Qv||- The condition of Theorem 12 . II admits the 
following weakening. 

Theorem 3.1. The process £ (t) is stable if and only if 



lim 

N— >oo 



KAN) 

N 



lim — ||Qjv| 







Proof. The Hilbert norm \\A\\ of any square matrix (or any nuclear operator A) 
satisfies the inequality 

\A\L 



(11) 



where WA^ and \\A\\ 2 are the nuclear and Hilbert-Schmidt norms, respectively. 
Recall that || A\\ for positive definite A coincides with the maximal eigenvalue X m (A) 
of A and for A > {llj takes the form 

TrA 2 

(12) K (A) > w . 

The application of (|12|) and ijHJl to Qn by virtue of the equality TtQn — Nb (0) 
gives 



N 



(13) 

On the other hand, 

Hence, 
(14) 



lim^> 

N ~ 2 



K (N) < \ TrQ 



'N 



Ihn^)< 

W->oo N 




□ 



The inequalities (fT3"|l . fbl| can be specified. 

Theorem 3.2. Given the sequence of maximal eigenvalues (norms) {X m (N)} of 
positive definite Toplitz matrices {Qn}, generated by a non-negative measure da (9) 
10), it holds that 

K (N) 



lim 

N— too 



N 



maxm„ 



Proof. The Toplitz matrix Q n generated by the non-decreasing function is the 
sum of non-negative Toplitz matrices Q$ and QfP , generated by non-decreasing 
functions a c and ad, respectively. Let us denote by Am (N) and Am (N) the 
maximal eigenvalues (norms) of the matrices Q$ and Qa? , respectively. Since Qn 



> Q 



(d) 

N ' 



then 



A£? (N) < A m (A) = 



o {d) 

Wn 



Q (c) 



n 



< 



Q (d) 



!N 



o (c) 
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By virtue of Theorem 13. II A^ (N) — o(N). Hence it remains to prove that 
15 nm = maxm„ . 

V ' N^oo N a 

To this end, let us consider first the case of ad (0) having only a finite number 2s 
of points of growth. We will not use in this proof the fact that the jump points 
of (Td (9) are located symmetrically with respect to the point 8 = 0. The Toplitz 
matrix = (pd (j — k))^2o > 2s <C N, generated by ad, can be represented in 
this case in the form 

2s 

(16) Q$ = y^m a (-,e a )e a , 

Q = l 

where 

(•,e Q ) e a = (exp (i (j - k) a ))f^=o 
are N x N matrices of unit rank, so that Qn transforms a N x 1 column vector 
x = ( x i)^=v int0 

2s 

(17) Qn* = X! m « ( x ' e ") e " ' 

a=l 

where (•, •) is the scalar product in the linear space of N x 1 column vectors Cn 
defined in a standard way: 

JV-l 

(x,y) = x jVv x = fe)^) 1 > y = iVi^o ■ 

3=0 

Notice that the vectors {e Q } are linearly independent. Indeed, suppose that 
there is a set of complex numbers {z a } such that 



(18) J2 z <* e <* = ■ 

a=l 

Due to l|18[). the numbers z a satisfy the homogeneous system 

s 

exp (ik9 a ) z a = 0, k = 0, 1, 2s — 1. 

a=l 

But the determinant of this system is the Van der Monde determinant, which 
vanishes if and only if among the numbers {exp (ik6 a )} there are equals. The 
latter is impossible by our assumption. Hence all z a = 0. 

Let A be a non-zero eigenvalue of and hA be a corresponding non-zero 
eigenvector: 

s 

(19) Y m a (h A , e„) e a = Ah A . 

a=l 

By (|T§|| h\ admits the representation: 

s 
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where z a are some complex numbers, not all of which are equal to zero. Put 

Va = \/"^( h A, e a) • 

By virtue of (|19l) . not all numbers r\ a = 0. Taking the scalar products of both sides 
of (|19|l with all vectors ^/m Q e Q , we obtain the following homogeneous system for 



(20) 



a'=l 



Thus, the non-zero eigenvalues of Q$ coincide, with account of their multiplicities, 
with the eigenvalues of the 2s x 2s Hermitian positive definite matrix 



(21) A N = {^Jm a m a > (e a ,, e a ))£ a , =1 

Notice that by definition of the vectors e a , we have 



(22) 



exp (iN (9 a > 



0)-i 



, a ^ a; (e a ,e a ) = N 



exp (i {9 a i - 6 a )) - 1 
Hence, the matrix An is the sum of the diagonal matrix 

4i,jv := (Nm a 6 aa <) 2 a 3 a , =1 

and the Hermitian matrix A2.N with zero diagonal elements and non-diagonal ele- 
ments y/m a m a i (e a ',e iQ .) , a'. By (1221 the non-diagonal elements of A2,n are 
uniformly bounded: 



\y/m a m a > (e Q /,e Q )| < 2 ^max|0 Q / — 9 a \ 1 ) (maxm Q ) , 



and, hence, 



\A 3 ,n\\ < (4s 



2) ( max \ 9 a ' — 9 a \ 1 ) [maxm Q 



Therefore, 



N - (4s - 2) max \6 a > - 9 a 



\A 



2,N\ 



< 



maxm„ < \\A 1<N \ 



iV+(4s-2) (max\9 a , -9 a \ 



I -4 



2, AM 



< 



( maxm a 



H-AjvII < \\Ai, N 
We see that 

(23) \$(N) = \\A N \\ = TV • maxm Q + (1) . 

N—too a 

To prove the relation (|15|l for a non-decreasing step function <7<j (6*) , (71-) — 
( — tt) < 1, having infinitely many points of jump, we take a small £ > and 
split ad (9) into a sum aid {9) + a^d (6) of non-decreasing step functions a\ : d (9) and 
&2,d (9), where, as before, aid (9) has a finite number of jump points and a2d (9) is 
such that 



/da2 d(9) < £ < maxm a . 
a 
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With respect to this split, we represent the Toplitz matrix Q$ as the sum Q^'^ + 
Q^' d ^ of non-negative Toplitz matrices generated by (9) and o^d {&)■> respec- 
tively. Notice that by construction 

(24) 

Besides, 



n (2,d) 



< TrQ%' d) < Ne . 



IN 



Applying the estimate lll'-'il) to and taking into account the inequality 124|) 

for N — > oo yields 



N ■ maxm a + 0(1) = \£ d) (N) < \$ (N) < N ■ ( 

a 

Finally, 



maxm Q + e 

\ a 



O(l) 



, ,. AL d) (N) . — \Z> (N) 
maxm„ < iim < lim < maxm Q + e 

N ~ N^oc N - a 



N- 



where e > can be taken arbitrarily small. 



□ 



Remark 3.3. For the Hilbert norm ||C|| of any N x N matrix C — {cjk)^ k l a the 
following estimate holds: 



N-l 



N-l 



(25) ||C|| < max < max \ c jp\ i max \epk\ } ■ 

{ 3 P =o fc p=0 J 

If C = [cj-k)j^l is a Toplitz matrix, then by l|25(l 

N-l 

(26) ||C|| < | Co | + J2 (M + Ic-pI) ■ 

P =i 

Thus for the Toplitz matrices Qn, which are under consideration here, Theorem 
I3.2l bv virtue of iffify . gives 

1 1 1 W_1 

(27) -maxm Q = lim — ||Qjv|| < lim — V \b(p)\ . 

Therefore the condition 
(28) 



p=0 



N-l 



is sufficient for the process stability. 



Jim l£>(p)|=0 



p=0 



4. Extension to continuous time processes 

Real signals are, certainly, continuous time processes, £ (t). The correlation 
function b(t) of a process having a finite second moment (£ 2 (t)) is a Hermitian 
positive function. As such, it admits the representation 

/oo 
exp (iAi) di? (A) , 
- OO 
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where i? (A) is a bounded non-decreasing function on the real axis. Like for the 
discrete time processes, i? (A) can be represented, in general, as the sum 

$ W = fie (A) + dd (A) 

of a non-decreasing continuous function # c (A) and a non-decreasing step function 
§d (A), and we call the process stable if $ (A) is continuous and unstable otherwise. 
To investigate the instability characteristics of a continuous time process, we con- 
sider instead of the Toplitz matrices Qni the set of non- negative integral operators 



(30) 



(B T f)(t)= f b(t-s)f(s) 
Jo 



ds, < T < oo , 



in the Hilbcrt spaces L 2 (0, T). Since b(t) is a continuous function, all these oper- 
ators are nuclear and their nuclear and Hilbert-Schmidt norms ||-Bt|Ii an d ||-By|| 2 
are given by the expressions 



(31) 



\B 



Tib 



Tb (0) 



<W(A) 



Bt\\ 2 = 



\ 



2Tjf \b{t)\ 2 dt = 



(32) 



\ 



(A - A')' 



■ sin 2 ^ P T dti (A') di9 (A) 



Let us denote, as before, by {m a } the set of jumps of i?(A). Using (f5T|) . i|3"2jl and 
the arguments similar to those employed in the proofs of Theorems 12.11 and 13.21 we 
obtain the following criteria of stability of a continuous time process £ (f ) . 

Theorem 4.1. A stationary continuous time process £ (f ) is stable if and only if 
its correlation function b (t) possesses the property: 

T 



Otherwise, 



lim — 



T 



Theorem 4.2. A stationary continuous time process £ (i) is stable if and only if 
the Hilbert norms \\Bt\\ of integral operators where 6(f) is the correlation 

function of the process, are such that 



lim - \\B T \\ = 



Otherwise, 
(33) 



maxm 
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Remark 4.3. For the norm of the integral operator Bt the following estimate: 

T 

||Br||<2 J \b(t)\dt 
o 

holds. As it stems from (|33|) the relation 

T 

lira - ( \b(t)\dt = , 

T^ooT J 


guarantees the stability of the process £ (t) . 

5. Application to processing of random signals 

5.1. Detection of quasi-periodic components in a random signal. The jumps 
of the spectral distribution function i? (A) of a stationary process is, in general, a 
sign of appearance of undamped oscillation components in a signal and the points 
of discontinuity of (A) are either the frequencies of such components themselves 
or directly related to them. Notice that, due to physical reasons, the measurement 
of £ (t) is possible only at discrete moments of time with a step A. If t in (|29|l is 
an integer multiple of A, then it is evident that 

(34) b(t) = / exp (M) da {9) , ft = ^ , 

J-n A 

oo 

a (9) = [# (0 + 2nfi) - i? (2nft - ft)] , -ft < (9 < ft . 

n— — oo 

The function a (9) is bounded and non-decreasing in the interval [—ft, ft]. If t? (A) 
loses its continuity at the points ±Ai,±A2, ■ ■ • , then a (9) has a non-void set of 
discontinuity points 

(35) H=(E(t + J)4)2fi}c[-fi ) 0], 

where E (x) is the fractional part of the number x. (In general, ±#ji coincides 
with every ±^ 2 such that A^ — Xj 2 is a multiple of 2ft.) Therefore, in general, the 
jump of a (9) at a point 6'j is the sum of the jumps of # (A) at all co-images of flo- 
under the mapping (|35|1 .) Taking A as the time measurement unit, we return to the 
representation of b (t) for integer t. Thus, the spectral distribution function u (9) 
inherits all discontinuities of d (A) from the interval [—ft, ft] and also may get new 
ones at the points (calculated according to (123) related to the discontinuity points 
of i? (A) outside this interval. We see that the spectral distribution function for the 
discrete time process obtained in such a way from a continuous time process, has 
a non-trivial component ad if and only if the corresponding spectral distribution 
function of the initial discrete time process has non-zero jumps on some set of 
points. In other words, the values of a random continuous time process measured 
at discrete moments of time, form a stable discrete time process if and only if the 
initial process is stable. 
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The correlation function of the discrete time process delivers not only the de- 
scribed gauge of instability of the process, but also the following tool for the detec- 
tion of the points {±6 a }, which are the discontinuity points of a (0). Put 

N-l , , v 

S N (0) = b (0) + 2 ( 1 - ) b ( k ) cos k9 = 

k=0 ^ ' 



, s f sin 2 iJV (9 - 9') 

36 = / - 2 \ —da (9) 

I iV 2 sin 2 ± (0-0') y ' 



It is not difficult to sec that 



(37) lira 4©w (0) = v (0 + 0) - a (9 - 0) . 

TV— »oo N 

Further, take a sufficiently large N ^S> 1 and split the interval [—n, ir] into equal 
segments of longitude S such that NS < 1. Let <r (#) have a jump m a within the 
interval (IS, (I + 1)6). Since 

|N < M (i-y) < l sina; l < M . 

then, for 16* — IS\ ~ (5, we have 

eAr (0)> / "T ** { - 9 ~ 9 \ da(&)>-N I da(0')>-Nm a . 
y ' ~ J N 2 sin 2 4 6» - 9') V ; " 6 J V ' ~ 6 

Z<5 ^ IS 

On the other hand, if the continuous part <j c (9) of a (9) satisfies the condition of 
Proposition 12 . 21 then one can use the arguments similar to those employed in the 
proof of this proposition to show that 



Q N (0) = O (n*&^ 



at the points remote from the jumps of a (0). Hence, the values of 0jv (IS), where 
I is an integer which satisfies the condition — ir < IS < tt, at the distances ~ S from 
the jump points of a (0), must be visible as larger than those at the distances^ loS, 
where NloS ~ 1. 

The assertion of Theorems 12.11 and 13.21 can be used for the detection of symp- 
toms of emerging instabilities of a random process, which can be considered as 
stationary for long time intervals. The method consists in the construction of the 
correlation function of the process from a piece of its time series from the beginning 
of observation to a rather far off moment of time T in the future. Set, as usually, 

T — k T 

b(k) = Y~k^t {p + kHip) ~ m2 > m = ' 

J9=0 p=0 



and compute, for a sufficiently large N <C T, the numbers 



JV-l / , \ ,2 /. n -. AT-1 



(38) ij HQ^a _ 1 + 2 £ fi-A^Lg) 



fc=l V / W p = Q 



Hp) 



b(0) 
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or the numbers 

T T 

iy(i-|) \b( t )\ 2 dt, i y \b(t)\dt 

' 

for a continuous time process. An explicit tendency of any of these numbers to be 
bounded, for increasing N, from below by certain positive numbers, is a serious 
evidence of the process instability. 

The following example demonstrates that the manifestation of such a tendency 
begins the sooner in N the larger is the contribution of the oscillating components 
generated by d<Jd {9) into b (0). 

Let the correlation function of a stationary random process be given by the 
expressions: 



(39) 6(0) = 1; b (k) = b (-k) = ^m a cos k6 a , k = 1,2,. 

a=l 

s 

1 < s < oo, < m a < 1, < Ox, S < 7T . 



a=l 



The spectral distribution function a (9) of such a process is the sum of the spectral 
distribution function a c (9) of the "white noise", 



da c (9) = — d6, p = 1 - y^m„ . 



2^ 

a=l 

and the step function ad (0), the jump points of which are {±6*0}, and 

o d {-9 a + 0)-a d {-9 a -Q)=tr d (9 a +Q)-tr d (9 a - 0) = ^m c 

To calculate the right-hand side of (JSJ) in this special case, notice that 

1 } sin 2 \N (9 - 9') , 

1 2 V ' d9' = N . 



2tt J sin 2 \{6- 9') 

— 7T 

Therefore we see that now &n (0) = pN~ l and thus 

a=l 



1 / sin 2 |iV (6> Q - 6 a ,) t sin 2 ±7V (6> Q + 6 a .) 



Hence, the first term on the right hand side of (|40|l becomes dominant for 
N>2p(2-p) (l>a) • 
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FIGURE 1. Sequences of HOatH^ f° r the signals c^JprmS and 
c4Jprm22, linear scale of N. 



5.2. Numerical results. Let us apply the latter results to the investigation of real 
signals obtained from the Forsmark 1&2 boiling water reactor (BWR) VGM01 , 
and also to simulated random signals. 

In Fig.l and Fig2 we display the sequences of the values 



(41) 



^2 WQnWI 



N-l 

N 2 51 

j,k=0 



for the Toplitz matrices constructed for two different shots of real signals obtained in 
the Forsmark BWR, named in |VGM0l) as c4-lprm3 and c4-lprm22. The sampling 
interval of these signals equal .08 s, and both of them consist of 4209 points. In 
VGM(H] several methods were used to obtain the corresponding decay ratio (DR) 
value. The decay ratio is a parameter related to the signal stability: the signal is 
more unstable if its DR is higher, see Appendix for details. The mean values of 
the DR given by these methods are .90 for c^JprmS and .51 for c4-lprm22. 

We observe the difference between a more unstable shot, c4-lprm3 and a more 
stable one, c4-lprm22. As expected, for higher DR, the sequence tends to zero 
slower. 

To determine the points of discontinuity of a (9) we have also calculated the 
function 



N 



On (0) 



HO) 
N 



N-l 

-Y 

k=0 



k 

N 



b (k) cos kO 
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Figure 2. Sequences of || <5jv || 2 f° r the signals c4-lprm3 and 
c4-lprm22, logarithmic scale of N. 



for different values of N. The corresponding results are provided at Fig. 3 and Fig. 4, 
they demonstrate different behavior of jjQn (0) with N growing, and for different 
degrees of stability measured by DR. In both cases the number of segments into 
which the interval [-tt , tt] was split, was equal to 3000. It can be observed that for 
c4-lprm22, a more stable signal, the function jj&n {&) tends to zero with increasing 
./V much more rapidly. 

The peaks of the signal correspond to the instability frequencies. The main peak 
for c4-lprm3 is obtained for 9 = 0.24 rad, which corresponds to the frequency 
of / = 0.48 Hz. On the other hand, for c4-lprm22 the main peak is located at 
9 = 0.27 rad (/ = 0.53 Hz). We can compare these results with those obtained in 
VGMOT]. The range of main instability frequencies obtained there is [0.480,0.495] 
Hz for c4-lprm3 and [0.450, 0.557] Hz for c4-lprm22. The method allows to locate 
secondary instability frequencies situated quite close. In Fig. 3 there is a secondary 
peak located at 9 — 0.27 rad, it coincides with the main peak of the shot c4-lprm22. 
Notice also that the present approach to the detection of instabilities is model-free. 

In addition, we have also simulated signals with different degrees of stability 
using the model suggested in jS VM88j . i.e., from the following continuous Langevin 
model: 



(42) 



at)+ci(t)+u(o=F(t), 
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with F(t) being a Gaussian colored external force such that 
(43) F{t) + T- 1 F(t)=T- 1 W(t), 

and 

U(0 =a^ + a 2 e + a 3 e. 
Here c is a damping constant, cij, j = 1,2,3, are some model constant parame- 
ters, and t is the correlation time, while W (t) is a white Gaussian noise with the 
correlation function 

(W(t)W(t')) = DS(t-t'). 

The parameters of the model, c, a\, 02, 0.3, D and r, are directly related to 
the stability of the signal: a\ coincides with u> 2 , w = 27r/ being the fundamental 
frequency, and the parameters c and a\ determine DR, see Appendix. 

In order to compare the results obtained for the simulated signals with those 
obtained previously for the real ones, we constructed the simulated signals with 
the values of the parameters chosen to produce the values of / and DR similar to 
those of the signals c4_lprm3 and c4-lprm22, and with the same sampling interval, 
0.08 s, and the same number of points, 4209. In Fig|S]and Fig[fJlthe sequences of 
the values of the Toplitz matrices for the models with c = 0.689 (DR = 0.5) and 
c = 0.105 (DR = 0.9) are gathered (the values of other parameters are: ai = 9 
.87, which corresponds to / — .5 Hz, a 2 — 03 = 0, D — 500, r = .6). 

As expected, the behavior of these signals is close to that obtained for the real 
signals, the sequences tend slower to zero for higher DR. 
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FIGURE 4. Function ^0jv (#), for the signal c4-lprm22, N = 100, 
N = 300, N = 500. 



For the same simulated signals, the discontinuity points are shown in Fig. 7 and 
Fig. 8. As before, the number of segments in which the interval [-tt , tt] was split 
was equal to 3000. 

Again, for lower DR the function jjQn {&) tends to zero with N much more 
rapidly. The main peak is now obtained at 6 = .24 (/ = .48 Hz) for the model 
with DR = .5, and at 6 = .25 rad (/ = .50 Hz) for the model with DR = 
.9. 

5.3. Appendix. The notion of decay ratio (DR), a basic parameter in the analysis 
of reactor stability, is deduced from the oscillatory model H42JI . but with STV97 



The DR is a measure of the system damping defined as the ratio between two 
consecutive maxima of the signal, for the model H42|) . 144|) it is a constant parameter, 



Neutronic signals are very noisy and, in general, their behavior cannot be fitted to 
that of a continuous second-order system, hence the DR in reality is not a constant, 
and its value depends on the model used to evaluate it VGM01 . Nevertheless, it 
gives a hint on the system stability. 



(44) 
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Figure 5. Sequences of ||<2jv|| 2 f° r the model signals, linear 
scale of N. 
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Figure 6. Sequences of j^j- ||Qjv||2 f° r tne model signals, loga- 
rithmical scale of N. 




Figure 7. Function -^Q N (9), for the model with DR = 0.5, 
N = 100, N = 200, N = 400. 
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Figure 8. Function j r Q N (9), for the model with DR = 0.9, 
N = 100, N = 200, N = 400. 



